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ABSTRACT AND SUMMARY 


Hanson (J. Acoust. Soc. Am., 54, 6, 1973) has modeled turbomachinery 
noise in a free-field as a sequence of pulsed-dipole acoustic sources with 
randomly modulated amplitudes and arrival times. This approach produced a 
far-field noise spectra with a broad-band background in addition to the 
usual discrete noise component. A systematic method for calculating the 
internal and external acoustical fields of aircraft turbofan engines has 
been developed by Zorumski (NASA TR R-419). This work shows that the duct 
acoustic problem can be reduced to a linear matrix equation WA - Q where 
A and Q are vectors of modal amplitudes and V is a square transmission/ 
reflection matrix. In this paper the methodologies of Hanson and Zorumski 
are combined to yield a unified analysis of ducted turbomachinery noise. 

The essential features of Hanson's analysis are unchanged by the addition 
of duct acoustic effects. It is shown that the far-field broad-band and 
discrete noise spectral components can be expressed in terms of modal cross- 
spectral matrices and directivity vectors which are derivable from the duct 


analysis. 



INTRODUCTION 


Aircraft turbomachinery noise has been analyzed by various authors 
with the aid of either of two possible simplifying assumptions. The noise 
originates inside the engine at the rotating arid fixed blades of the turbo- 
machinery stages, consequently it is possible to study the ducted sound 
field assuming that the engine duct is infinitely long. This was the main 
viewpoint used in the classic analysis of "spinning mode" generation by 
Tyler and Sofin (ref. l). For these authors, the radiation of the noise 
was considered as an independent problem. Using an alternate simplifying 
viewpoint, Lowson (ref. 2) argued that the main features of radiation 
from engine fans could be found by ignoring the duct effects entirely and 
considering the noise sources to be in the free-field. Both of these views 
are valid up to a point and may be used when only qualitative effects are 
needed. It should be recognized, however, that such drastic simplifications 
do not make sense in some situations. For example, if one is studying the 
effect of duct acoustic treatment on fan noise, it would hardly be sensible 
to neglect the duct; or, if it is desired to predict noise on the ground 
during an aircraft flyover, it is not sensible to ignore the radiation from 
the engine inlet. Recently Hanson (ref. 3) has used the no-duct assumption 
to great advantage to develop his "Unified Analysis of Fan Stator Noise." 
Hanson assumed that the essential features of his analysis would not be 
changed by the addition of duct effects and it will be shown in the present 
paper that his assumption is correct. In order to do this we will employ 
the results developed by Zorumski (ref. 4) for the radiated field due to 
ducted noise sources. Thus, the purpose of this paper is to combine the 
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analyses by Hanson and Zorumski to develop a unified analysis of ducted 
turbomachinery noise. 
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SYMBOLS 


a 


A 

b 

B 


c 


D (<£) , c (4>) 

E l } 

f (t) 


F (w) 


G 


I 

k 


m 


n 


P 

Q 


random variable with unit mean 
mode amplitude 

random variable with zero mean 
number of rotor blades 
ambient speed of sound 
mode directivity functions 
expected value of O 
force pulse 

Fourier transform of force pulse 

power spectral density 

intensity 
wave number, oj/c 

circumferential mode number 

harmonic number 

acoustic pressure 

source mode amplitude 
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r, 6, z 
r, 0, <j> 
B 


s 


t 

T 

V 

w 


X, y, Z 


Y 


Z 


6 (u) 



* 


cylindrical coordinates 
spherical coordinates 
mode reflection coefficient 


stator vane number 


time 

mode transmission coefficient, also rotor period 2tr/ ft 
number of stator vanes 


wave coefficient 


rectangular coordinates 
broadband noise directivity function 
harmonic noise directivity function 
Dirac delta function 


ambient density 
standard deviation 


radial mode function 
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to 


( ) 

[] 


radian frequency- 
rotor speed 
rov matrix 


column matrix 


square matrix 
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REVIEW OF HANSON'S ANALYSIS 


Hanson (ref. 3) considered the noise radiated from an aircraft 
engine fan stage due to the viscous wakes from the upstream rotor. Each 
stator vane source was modeled as a simple dipole (concentrated force on 
the fluid) which "pulsed" as each viscous wake from the rotor passed over 
the vane. The stator noise source model was therefore a circular dipole 
array, such as in figure 1, in the free-field. A typical pulse from one 
of these dipoles is f(t) and the transform of this pulse is F(oo). The 
pressure spectrum due to a single typical pulse at time t=o on stator 
vane s is 


p(r 



kr ■F(co) ikr 
s e s 

irrrr 2 

s 


( 1 ) 


Hanson recognized that the viscous wakes from the rotor blades vary in a 

random fashion and modeled the actual pulses on the stator vanes as having 

random amplitudes, with unit mean and standard deviation ^ a . The time 

of the pulse occurrence was modeled as a Gaussian process with zero mean and 

standard deviation . Figure 2 shows a typical sequence of pulses with 

this random amplitude and arrival time property. With this model, the 

intensity per unit bandwidth dl at the observer position 0 in the far 

du> 


field was found to be 
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dl 

dco 


O 14 I p 

Bfloo I F (h))f 
8 n 2 £> c 3 r 2 


C 


i 


<Ta + 1 - 


-(a® or 


> 2 ] 


6 ) 


broadband term 


+ Bft 





(oj-nBfi) E (<f>,0;u)) 


( 2 ) 


harmonic terms 


In equation (2), B is the number of rotor blades, V is the number of 
stator vanes, ft is the rotor angular speed, and T = 2lT/ft. The term 
containing GT in equation (2) is the broadband noise due to pulse ampli- 

3 , 

tude modulation (PAM)and the term remaining 1 - e is the broadband 

noise due to pulse position modulation (PPM). Note that PAM generates broadband 
noise without changing the harmonic noise, but that PPM decreases the har- 


monic noise by the factor , 

In the following sections, we will develop a formula, analogous to 
equation (i) ,for the far-field pressure due to a dipole pulse in a finite 
length duct and then use this result to extend Hanson's analysis to include 
duct effects. It will be shown that the essential results of the unified 
analysis are unchanged by the addition of these duct effects. 
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DIPOLE IN INFINITE DUCT 


In order to develop a formula for the far-field pressure due 
to a dipole pulse, we first consider the sound field due to a dipole in 
an infinite duct. In this paper, only hueristic arguments will be 
presented. The analysis of reference U is being followed here, so that 
the reader who is interested in analytical details may refer to that 
work. 

Consider the infinite duct which is formed by an imaginary ex- 
tension of the finite duct as shown in figure 3. The planes z and z 

12 

denote ends of a duct section, and a single dipole is located at 0=0, 
r-r in the source plane z . In axisymmetric ducts, it is always con- 
venient to develop the pressure spectrum in a Fourier series 


pc 

P (r,6,z;w)= — 



im0 

e 


p (r,z,w) 
m 


(3) 


and to treat each term in this series independently. Without bothering 
with a lengthy derivation, it is assumed that the sound field at z^ is 
made up of circular duct modes traveling away from the source to the 
right. These modes, which are denoted by ^ (r) are known to be pro- 

portioned to Bessel functions of the first kind of order m. The subscript 
m is the radial mode index. At z 1 then, the m til pressure harmonic in 
equation (3) is 
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p m (r, V c) - * 3 J2 . t^ - L jP % (kr) (4) 

( m»l 

In equation (4) is the source mode amplitudes at z and the positive 

superscript is used to denote the positive direction of travel. In practice 
the infinite series in equation (4) is alvays truncated so that it may he 
replaced by the matrix formula 


P m (r, V w) - L* mu (kr)J j<£jJ 



(5) 


The inner product of the row matrix and the column matrix in equation (5) 
replaces the summation in equation (4). 

Equation ( 5 ) gives the pressure at z ± due to a dipole if the duct 
is infinite. After considering some elementary transmission and reflection 
coucepos, we will show that this information can be used to predict the 
pressure in a finite duct. 
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TRANSMISSION AND REFLECTION 

It is known that the propogation of waves in uniform ducts is 
described by their axial wave number through the function e ik zmp Z ; 

that is, the amplitude of a mode at z is related to the amplitude of 
that same mode (in a uniform duct) at z by the equation 


A +2 

my 


ik 

e zmy 




( 6 ) 


In general k is a complex number so that equation (6) describes the 

phase shift and the attenuation of the wave between planes z and z . 

1 2 

The attenuation, an exponential decay of the amplitude, is depicted 

graphically in figure k which shows the relative amplitudes of the waves 
at planes 1 and 2. All transmission effects in ducts, both uniform and 
nonuniform, are described by the matrix equation 


a +2 I = I~t +2+1 
my myv 

L - 


c 


( 7 ) 


In the simple case of a uniform duct, the transmission matrix 
a diagonal matrix whose non-zero elements are given by the exponent 
factors in equation (6). In nonuniform ducts, a more general expression 


M is 

L ~ mvv I 

mential 
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must be derived for the transmission matrix, however its basic form will 
be as shown here. Of course, a similar equation is possible for waves 
traveling from right to left in the duct. 

When acoustic modes travel to the end of a finite duct, there 
is partial radiation of the sound to the far-field and partial reflection 
of the sound as waves traveling backward from the duct opening. This 
reflection process is shown graphically in figure 5 where a positive 

*f 2 _ Q 

moving wave is reflected as a set of negatively moving waves 

The analysis of this radiation and reflection process is complex and cannot 
he treated here, A general solution for radiation from a flanged annular 
duct is given in reference 5 Bud the problem of radiation from an unflanged 
duct is treated in reference 6. Here, we only note that the solution to the 
radiation problem also yields the reflection equation which must take the 
form 


( & ~2 
)A . 


-2+2 

R 

V * 2 

j mu 


mpv 

| mv 


(8) 


Considering the simultaneous effects of sources, transmission, and reflections 
gives the following equation 


I 0 
-T I 

' -R 0 

\ 0 0 


A +1 

J?x\ 

= k 3 lT (w)l ' 

lO 

0 

V o 

0-R 

1 ~I-T 
! 0 * 

i 

A J 

A 1 

i 

pc ! 

o'" 
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Equation (9) will "be denoted by 



( 10 ) 


therefore, the solution for the mode amplitudes is 



= k 3 1 Lad 1 

r c 



-i 



(ii) 
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RADIATION FROM DUCTED DIPOLES 


The solution for the sound field inside of the duct has been 
developed in terms of the mode amplitudes at the duct ends. These mode 
amplitudes also give the radiated acoustic field as shown in reference 5 
and 6. Here, each mode which strikes the duct end is assumed to have a 
known radiation directivity function D (<J> 2 ) which has been determined 
from the solution to the radiation problem. Thus, the far-field pressure 
at the observer point 0 due to a single mode at the end of the duct is 


P m (r, <|>;(*)) 


2 

pc 

OJ 




(<P 2 ) 



( 12 ) 


and the pressure due to all modes at the end of the duct can he found, by 
superposition, to be 


p(r ,e,4>,to) = 




Substituting the solution for the mode amplitudes, equation (11) into 
equation (13) gives 


p(r,e,<J>,co) = kr lt( M ) l 


e itr 2 



im0 


(lk) 


m=-°° 


lh 



where 


L A 


D 

m 


(<fr 2 )J = L 


V* 2 >- 



( 15 ) 


Equation (l4), the far-field pressure spectrum due to a dipole in a duct, 
replaces equation (l) when there is an array of dipoles in a duct with ran- 
dom amplitude and pulsing times- Although this expression is more complex 
than the free-field formula, Hansons analysis can still he canned through 
with no essential changes. 
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RESULTS OF UNIFIED ANALYSIS 


Essentially, all of Hanson’s results are preserved except that the 
expressions for the directivity effects are more complex. In the ex- 
pression for the far-field broadband intensity, Hanson's directivity function 
is 


Y (<J>) 


sin^B cos^<{) +{l/2)cos^3 sin“h£ 


(16) 


provided there are more than two stator vanes. The angle 3 in equation 
(l6) is the mean inflow angle to the stator which determines the orientation 
of the dipole with respect to the plane of the stator vanes. The analysis 

with duct effects shows that the broadband directivity for the ducted 
stator is 


oo 


Y<0; 0 )) = (4tt)‘ 


E 

m= 


|D <4»|[|Q m i |Q? 

► i1(d* (0)1 

>- mv J h mp j L r 

nvJJ y mu | 


( 17 ) 


that is, the broadband directivity is axisymmetric but depends on the 
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directivities of the cylindrical mode orders m and radial mode orders y 
v as well as the modal source amplitudes ^ which the unit dipole would 

generate in an infinite duct. The square matrix in equation (17) is 
the modal cross spectrum for a unit dipole in an infinite duct. Attempts 
have been made by Bolleter and Crocker (ref. 7 ) and by Harel and Perulli 
(ref. 8} to measure this cross spectrum for general, sources. These 
measurement attempts have been partially successful. The present analysis 
therefore provides a link between the results of nonreflecting duct 
measurements » which can be made in the laboratory, and the prediction 
of far-field radiated noise with a finite length duct. 

Hanson also found a fairly simple expression for the directivity 
of the harmonic noise from a free-field dipole array. This was 


1 XT' 2 1 

V 

1 /JYis (e>e ^» v) + 

x /Yps ^* e » < ( , » v ) 

( S=1 j 

s=r 


where 


K 1 

r / \-i 

sin ip j 

'h' 

i i 

= sin8 cos <p - cos 0 sin $ sin 10- 

s [ 

| L \ V /J 

COS iJj. I 

si 


( 19 ) 


or 
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^s " 


nBM sin <j> cos f 0 


( e 


2tts 

V 



2imB | , 


( 20 ) 


and M is the rotational Mach number at the dipole radius. Note that the 
directivity of the harmonic noise is not axisymmetric. The directivity 
of the harmonic noise from the ducted dipole array depends on the properties 
of the Tyler-Sofrin spinning modes. It is shown in the appendix that 
this expression is 


OO CO 

zi<t>, e ; w) = (47r) 2 X) £ 

k” k = -oo 




mp ) 


J(k'“k)V 0 

* 

( 21 ) 


where m = nB - kV, (22a) 

and m 1 = nB - k f V ( 22b ) 


The harmonic directivity depends on only the selected mode orders m=nB-kV 
as predicted by Tyler and Sofrin; however it is not axisymmetric because of 
the factor exp [i (k -k)V6 j which corresponds to Hanson’s result for the 
free-field array. If the intensity is averaged over the azimuth angle 6, 
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the double sum in the directivity expression (21) is reduced to a 
single. sum. 
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• CONCLUDING REMARKS 


A unified analysis Las been formally carried out for the radiated 
noise field of a circular array of pulsing dipoles in a finite length 
duct. The result shows that there is no qualitative difference in the 
predicted broadband and harmonic spectra between a free-field array 
and a ducted arrays however, the magnitude and directivity of the radiated 
field is changed by the presence of the duct and a quantitative computational 
approach has been outlined which will predict the duct effects. The 
analysis has also shown that measured modal cross spectra from sources in 
a nonreflecting duct may be used to predict the noise radiated from those 
same sources in a finite duct. The analysis also shows again that the 
harmonic noise is made up of the familiar Tyler-Sofrin "spinning modes" of 
order m=n£-kv . 
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APPENDIX 


Statistical Analysis 

In the text a relation is presented (eqs. (lit) and (15)) for the 
far-field pressure spectrum (complex pressure) due to a dipole pulse at 
position = 0 in an acoustically lined duct (fig. 6). If the pulse 

occurs at time t = 0, the far-field pressure is 

fi(r,6>^) - (xr) 1- %^ — jr LD^(0)j[.w rri i l {Q w } ^(e-e s ) (A-l) 


However, if the pulse occurs at an arbitrary time t = t the pressure 
transform is 


Lcy£ 

e 


o 




Now consider the complex pressure p (r ,0 ,<f> , to) produced in the 

s n 

far-field due to the rr^- pulse on the s— vane. The amplitude of 

this pulse will be written as a |f(w)| where a is a random 

sn sn 

variable with unit mean, that is 


E n] - I 


(A-2) 


for all values of s and n. Therefore, as discussed by Hanson (ref. 3), 

J F(to) | is the mean pulse shape and the random variable a accounts 

sn 

for the pulse amplitude modulations . The reader is cautioned that Hanson 
uses the subscript m in reference 3 as the vane number, whereas, in the 
present paper s denotes vane number and m denotes the circumferential 
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for the n— 


vave number in the duct analysis . The arrival time 
pulse on the s vane, as developed by Hanson, is 


t 

sn 


*sn = ( Sl - & 

W B 



4 - 



T 


(A-3) 


where V is the number of vanes , B is the number of blades , T = — ,• c 

ft 15 

the period of rotation s ft being the rotor angular velocity. The parameter 
b sn iS a rand - om variable which accounts for the jitter or pulse position 
modulations. Here it is assumed that 


e {>«*} = o 


(A-4) 


for all values of s and n, 


Hence, the mean arrival time for the n^- pulse on the s— vane is 


- h. + (A-5) 

The parameter q s in equations (A-3) and A-5) is defined by 

(a-6) 

where the bracket [ ] indicates that only the integer part of the 

quotient is retained. Substituting equation (A-3) into equation (A-l) gives 



<fsn(r,e, 2 r,oj) - o. sn e ^ wTb sn T (Kr )lffc J )\ 




o(0-e Sj 0) 


(A-7) 


D(e>) = I Liu»)JLw m f ’{$„] e 

^ V -Pa 


.CKT© 


(A-8) 


An observer located at a fixed position in the far field will 

experience a pressure time-history produced by a sequence of pulses as 

shown in figure 2. Thus, a particular sequence of realizations for the 

random variables a and b determine a single realization of the 

s n sn 

random far-field pressure. Equation (A-7) gives the Fourier transform of 

a typical far-field pulse p (r,0,<f>,w) which itself is a random variable. 

s n 

In the following it will be shown that the power spectral density G(r 

of the far-field pressure can be determined in terms of statistical averages 

performed on the random variables p (r,0,<f>,w). An expression is then 

sn 

determined for the acoustical intensity per unit bandwidth. 
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First rewrite equation (A-7) as 


ifsnUe^w) - X iYi + l Ysn 


- 0-SY\ D5I 


e 


n.wCcsn - *' b sr \T ) 


(A-9) 


in which 


d s , irxiiSyiL D (©-©=, 0 > 

An = <27 

^■SY\ — A 4' ^5 


> 

-> 


(A-10) 
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The power spectral density for the far-field pressure p(r,0 5 <J>,t) » 
consisting of an infinite train of pulses a is given by (ref. 3) 




L ' ^ B 
^ (2n+'Jt 




(A-ll) 


where : 


X X sn+ W 


M 'V* 


which can immediately be rewritten as follows : 


* Mtr<a ' # ' w) - ^ E [ ,* 9 7 } 


zn ^ 




(A-12) 
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Herein it will be assumed, that the 2V(2II + l) random variables 

a sn an( i kg n a re independent. It then follows that the 2V(2N + l) 

random variables X and Y are uncorrelated. An additional 

sn sn 

consequence of the assumption of independence is that 

(^sn, ^sr>) = ( k.sn!> ^sn( ^sn ) (A-13) 

where ^ sn ^ a sn ’ se ^ 2V(2N + l) joint probability density 

functions for the a 1 s and b *s and f (a ) and g (b ) are the 

1311 sn sn sn sn sn 

respective individual densities. It is further assumed that these 
densities are independent of s and n (vane number and pulse number). 
Therefore > the a sn ,s an, i ^ sn * s are assoc i a t e fi with two classes of 
random variables a and b, respectively. The required probability 
density functions will be assumed as 


f srv (dsr) - ? 


] 


(A-Ht) 


28 



for all n and s. It is not necessary to specify a density for the 
random variable a. However, from equation (A-2), E{a} = 1 and its 
standard deviation will be noted at A Gaussian density is 

assumed for the random variable b. As indicated by equations (A-U) 
and (A-l4) , E{b} = 0 and its standard deviation is ^7^. 

Continuing the development from eqaution (A-12) the following 
veil-known result will be useful: 



(A-15) 


where VAR{ } denotes the variance of the expression within the 

brackets. Since the X t are uncorrelated 

sn s 


VftR 




(A-16) 
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where : 


Relations similar to equations (A- 15 ) - (A-17) can also be written for 

the Y variables . 
sn 

Substituting these results (eqs. (A-15) to (A-17)) into equation 
(A-12) gives: 




Z-TI 


r is/ V 

b)Z- X X s „ 

L p. * - w s -• 1 



-V 


t* 




i i h * 

ftS-N S* 1 


H V , w 

+ 1 II HY,*.} - l 


2£ 





H \r f N 2, 

- Y— *51 ( el { V 5K) ^ | 


(A- 18 ) 
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With the use of equations (A-13) and (A-lfc) the various 

expected values appearing in equations (A-18) can be easily calculated 
( ref . 3 ) and are 


= |Psl os(wc: stl )e 
e{Y S v,j - |d s \ e * 

£ ^ D s ^ £ L a ' 2 '} 




- Co T 






) (A-19) 


Substituting equations (A-19) into equation (A-18) gives 


- 2-n 1 




-2 2^ 
— c*j T Vb 


N V 
n * -M 


iD*\ + 


-WTV; r V 

^-zTie | { i- \o*\ cos\wc sy ^V •+ 

w -\r 

S Y\ 



(A- 20) 
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In equation (A-20) use has been made of the relation (since the random 
variable a has unit mean) 




(A-21) 


The power spectral density G(r ,0 ,<)>,w) (eq. (A-ll)) will be 
determined by separately evaluating the two terms in equation (A-20), 
Thus, let 
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It will be shown that G^(r,0 ,4 » ,oj) is a broadband component 
G 2 (r,9,<f>,to)is a harmonic (discrete) component. 

The double sum in equation (A-23) can be simplified to 
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whereas 
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and the following relation has been applied: 
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Now by substituting equation (A-25) into equation (A-23) and carrying 
out the limit as N -*■ « gives 
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which is the broadband component of the far-field spectral density. 

In order to determine the discrete component G^(r,d ,<f> ,<o) note 
that (ref. 3) 
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Now fcy using the results that (ref. 3) 
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the following result is obtained: 
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and 
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It should he observed that G (r^Q^suO is a discrete spectrum 
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can nov be expressed as 
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where 
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Figure 1.- Circular dipole array. 
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Figure 2,- Dipole pulse amplitude and position modulation. 






